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Salamanders of the genus Speleomantes are the only mem- 
bers of the family Plethodontidae occurring in Europe: sev- 
en species (S. ambrosii, S. italicus, S. flavus, S. supramontis, 
S. imperialis, S. genei, S. sarrabusensis) are endemic to Italy, 
while the distribution of S. strinatii also extends to France 
(LANZA et al. 2006). A few records of Speleomantes exist 
from beyond their natural range. Individuals of S. ambro- 
sii and S. italicus have been introduced to a cave in south- 
ern Tuscany (Italy) for an experiment of artificial syntopy 
(ForTI et al. 2005), which has resulted in the formation of 
new hybrids (CIMMARUTA et al. 2013). In other cases, in- 
dividuals of S. ambrosii have been released in the eastern 
part of the range of S. strinatii (Liguria, Italy), while indi- 
viduals of this latter species have been moved to the Italian 
Karst (northeastern Italy) (LANza et al. 2006). Outside It- 
aly, allochthonous populations of Speleomantes are known 
from Germany and France. In the first case, a population of 
about twenty S. italicus has been known from a quarry in 
a forested area in central Germany since 2013 (GINAL et al. 
2021, SCHULZ et al. 2021). One of the alien French popula- 
tions was discovered in central France in 2015, when about 
thirty individuals of S. strinatii were observed inside a nat- 
ural cave (LUCENTE et al. 2016). The second French alloch- 
thonous population was introduced to a mine in the Pyr- 
enees in 1970 (LUNGHI et al. 2018) in the shape of individu- 
als from a rearing program conducted at the Subterranean 
Laboratory of Moulis (Ariège, France) (DURAND 1967, b). 
The old registers of this experimental program report that 
the Speleomantes reared in Moulis had been collected from 


about fifteen localities ranging from the Maritime Alps to 
the Tuscan Apennines, an area that embraces the distribu- 
tion of all the three mainland Speleomantes species (GuIL- 
LAUME & DURAND 2003, LUNGHI et al. 2018). 

Considering the interfertility that characterizes the 
mainland Speleomantes spp. (CIMMARUTA et al. 2013, FICE- 
TOLA et al. 2019), it was recently hypothesized that the in- 
troduced population of the Pyrenees may consist of a new 
generation of hybrids (LuNGHI et al. 2018). This hypoth- 
esis was also supported by the peculiar dorsal colour pat- 
terns observed in some individuals, which seemed a mix of 
those characterizing the three mainland species (LUNGHI 
et al. 2020b). The scant existing information on the origin 
of the Pyrenean allochthonous population is mostly an- 
ecdotal and refers to S. strinatii (RAFFAELLI 2007); how- 
ever, no studies have been performed yet to shed light on 
the identity of this introduced population. In this study we 
assess the genetic structure of the Pyrenean Speleomantes 
population and compare it with the available information 
in order to unravel the identity of its components and ver- 
ify whether or not these cave salamanders represent hy- 
brids. 

In August 2018, we surveyed the allochthonous Speleo- 
mantes population living in a mine situated within the 
Parc Naturel Régional des Pyrénées Ariégeoises, in France 
(LuNGHI et al. 2018). As in previous papers, no further in- 
formation on the exact location of this site is provided to 
ensure the protection of this population (LUNGHI et al. 
2019). We collected tissue samples from 31 specimens of 


© 2022 Deutsche Gesellschaft für Herpetologie und Terrarienkunde e.V. (DGHT), Germany 


Open access at https://www.salamandra-journal.com 


67 


Correspondence 


Speleomantes captured inside the mine. Using a scalpel, we 
removed about 10 mm of the tail tip, placed it in a 2-ml 
plastic tube filled with absolute ETOH right after the cut. 
Before each cut, the scalpel was held over a flame until it 
became red hot; this ensured thorough sterilization and an 
easy cut. Total DNA was isolated from the tissue sample 
according to a slightly modified CTAB protocol (DOYLE 
& DOYLE 1987). Two mitochondrial regions of the genes 
cytochrome-b (cytb) and NADH dehydrogenase subunit 2 
(ND2) were amplified using the primers and PCR condi- 
tions as reported in CIMMARUTA et al. (2015). PCR prod- 
ucts were submitted to Macrogen Inc. (www.macrogen. 
com) for purification and sequencing. Electropherograms 
were visually checked and the obtained sequences were 
deposited in GenBank (accession number MZ666133). A 
524-bp fragment of the nuclear gene sodium-calcium ex- 
changer (NCX1) was also sequenced using the primers 
NCXı_ F (5’-GACTGTCTCCAACTTGACCTTGAT-3’) 
and NCX1_R (5-CTTTGAGGATTCTGGCCATGT-3), 
as originally reported in ROELANTS et al. (2007). PCR re- 
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actions were performed according to RoviTo (2010), in a 
total volume of 25 ul, containing 1 unit of Taq DNApoly- 
merase (Promega), 1 unit of reaction buffer 10X, 2.5 ul of 
MgCl2, 0.2 ul of dNTPs, 0.2 of each primer, and 2 ul of 
DNA template. Thermal cycling conditions were the fol- 
lowing: initial denaturation at 94°C for 5 minutes, followed 
by 35 cycles at 94°C for 30 s, 1-minute annealing at 51°C, 
1-minute elongation at 72°C, and final elongation for 10 
minutes at 72°C. The obtained sequences were deposited 
in GenBank (accession numbers MZ666135-MZ666137). 
The analyzed NCXi fragment was used to confirm the mi- 
tochondrial taxonomic assignment of the samples since it 
showed two diagnostic SNPs (positions 3 and 192) between 
the three mainland species of Speleomantes: TC S. strinatii, 
TT S. ambrosii, CT S. italicus. To assign the obtained se- 
quences to a species of Speleomantes, concatenated mito- 
chondrial haplotypes were compared to those of the three 
mainland species: S. strinatii, recovered from CIMMARUTA 
et al. (2015); S. italicus and S. ambrosii ambrosii, recovered 
from LucENTE et al. (2016), with S. imperialis and S. supra- 
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Figure 1. (A) Phylogenetic relationships between S. italicus, S. ambrosii and S. strinatii from concatenated sequences of cytb and ND2. 
Neighbor Joining tree; bootstrap values above 85% are shown on each node. The Sardinian species S. supramontis and S. imperialis 
serve as outgroups. (B) Median-joining network including 194 sequences of S. strinatii from CIMMARUTA et al. (2015) and the 31 se- 
quences from the Pyrenees analyzed in this study (squared). Haplogroups A1-4 and B1-3 are named and coloured as in CIMMARUTA 
et al. (2015). (C) An individual from the introduced Pyrenean population (LUNGuI et al. 2020b). (D) An individual from the identified 
source population located in the southern part of Genoa Province (LUNGHI et al. 2021); the dorsal pattern is similar to that from the 


Pyrenean population (LUNGHI et al. 2020b). 
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montis for further comparison. After finding the best sub- 
stitution model (TN93+G), these sequences were used to 
build a Neighbor Joining (NJ) tree based on TAMURA & NEI 
(1993) distances with set gene partitions and one thousand 
bootstrap replicates. All analyses were run in the software 
MEGA X (Kumar et al. 2018). A median-joining network 
was built on the mt concatenated alignment, consisting 
of the whole S. strinatii mitochondrial dataset from CIm- 
MARUTA et al. (2015) and the concatenated sequences from 
the Pyrenees, using Network v. 10.0 (http://www.fluxus- 
engineering.com) under the Greedy FHP criterion (BAN- 
DELT et al. 1999, FOULDs et al. 1979). 

The fragments obtained from the 31 individuals from 
the Pyrenees were 638 bp in length for cytb and 675 bp for 
ND2, providing a final concatenated alignment of 1313 bp. 
A single haplotype was found for the cytb gene; at ND2, 
two haplotypes were recovered that differed in a single mu- 
tation, with the same haplotype found in all the sampled 
individuals except one. Both NJ and ML analyses carried 
out on concatenated sequences provided the same topol- 
ogy, with the Pyrenean haplotypes nesting within S. stri- 
natii (Fig. 1A), and clustering with the samples from the 
environs of Genoa (southeastern Liguria). The same result 
was obtained in the outcome of the median-joining net- 
work (Fig. 1B), where the haplotypes from the entire spe- 
cies range of S. strinatii were included to represent all the 
subclades found in this species. One of the two haplotypes 
from the Pyrenees was the same as the most frequently re- 
covered one from the caves in the area of Rapallo, Bargagli 
and Isoverde (haplotype hst6), while the second haplotype 
(hst34) was a single mutational step away. All of the 31 Pyr- 
enean individuals showed the same bases (TC) of the other 
S. strinatii samples at the two NCX1 SNP positions. 

The genetic data obtained clearly identify all individu- 
als examined from this Pyrenean locality as S. strinatii. The 
mitochondrial haplotypes are identical to the most com- 
mon one found in southerneast Liguria (haplotype hst6; 
CIMMARUTA et al. 2015). The source area of the introduced 
Pyrenean population is therefore most likely located in the 
southern parts of Genoa Province, where the same hap- 
lotype is largely dominant and even exclusive in the local 
gene pool in the cave salamanders of the caves of Rapallo, 
Bargagli and Isoverde. The only differentiated haplotype 
(hst34) diverges from hst6 by a single mutation at the ND2 
gene. The analysis of the nuclear NCX1 gene fragment also 
confirms that the sampled individuals belong to S. strinatii, 
without any sign of admixture. 

The different haplotypes characterizing the two French 
allochthonous populations so far known (that showed 
up in this study and the one described in LUCENTE et al. 
2016) highlight two independent introduction events, 
with founders coming from different clades of S. strinatii, 
whose gene pool is subdivided into two main groups that 
inhabit the western and eastern parts of its range, respec- 
tively (CIMMARUTA et al. 2015). The individuals sampled 
in the Pyrenees and analysed for this study have haplo- 
types belonging to populations of the eastern clade, while 
those from Poitou-Charente province had haplotypes 


nested within the western clade of S. strinatii (LUCENTE 
et al. 2016). Among S. strinatii from the source popula- 
tions (Rapallo, Bargagli and Isoverde), patterns that di- 
verge from the “classic” yellowish coloration of this spe- 
cies are sometimes found (LANZA et al. 1995, LUNGHI et 
al. 2021) in specimens that look like the more “reddish” 
ones observed in the Pyrenean population (LUNGHI et 
al. 2020a) (Figs 1C-D). However, such reddish colora- 
tion seems to be more frequent and have greater intensi- 
ty in the allochthonous Pyrenean population, potentially 
suggesting the occurrence of stochastic events that have 
boosted the spread of a particular phenotype. On the oth- 
er hand, additional influences (e.g., environmental condi- 
tions or diet) cannot be excluded; indeed, S. strinatii has 
a relatively wide distribution range across different envi- 
ronments (LANZA et al. 2006), which is a condition that 
likely promotes phenotypic plasticity (AULD et al. 2010, 
LUNGHI & ZHAO 2020). 

The data presented here negate the hypothesis of a hy- 
brid Pyrenean Speleomantes population at the investigated 
site. Actually, none of the known hybrid populations has 
S. strinatii as a parental species (CIMMARUTA et al. 2013, F1- 
CETOLA et al. 2019). Furthermore, a study analysing an area 
of parapatric distribution between S. strinatii and S. am- 
brosii, as well as the results of an experiment of artificial 
syntopy in a cave in the same area, suggested S. strinatii as 
the more successful competitor and able to relegate S. am- 
brosii to the less suitable environmental sectors of the cave 
(CIMMARUTA et al. 1999). Therefore, even if individuals of 
all three mainland Speleomantes species were originally re- 
leased in the Pyrenees, S. strinatii would have outcompeted 
S. ambrosii and S. italicus in terms of adaptability, compe- 
tition, or even by its higher tolerance of suboptimal mi- 
croclimatic conditions (FICETOLA et al. 2018) and would 
therefore have persisted as the “last survivor” of the origi- 
nally introduced multispecies assemblage. 
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